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In this paper, the nonlinear theory of plasma waves is extended to the plasmas that their 
equilibrium state are specified by the non-Maxwellian (here kappa) distribution. We believe 
that the extension is very important since most of the space and some of laboratory plasmas are 
not in the Maxwellian equilibrium. Although the linear theory of this issue has been known for 
£\j [ the decades but, to our knowledge, this is the first attempt in opening the gate to the nonlinear 

world of plasma waves with the non-Maxwellian equilibrium. As an example the ion-acoustic 
solitons are studied in this framework taking into account the electron trapping in the trough 
of longitudinal field. It is shown, as the most important result regarding to the non-Maxwillian 
equilibrium, that there is the possibility for the ion-acoustic solitons to move subsonically. The 
solitons velocity and their width are monotonically increasing functions of spectral index (k) 
^ , and approach to the their Maxwellian values as n — > oo. 

I. INTRODUCTION 

CO 

O 

In the natural space environment, e.g., planetary magnetospheres, astrophysical plasmas, and the solar wind, 
plasmas are generally observed to possess a non-Maxwellian high energy tail [1-5] . A useful distribution function 
to model such plasmas is the generalized Lorentzian (kappa) distribution. Important features of n distribution 
£L|| are that, first, at high velocities the distribution obeys an inverse power law, and that, second, for all velocities, 
in the limit as the spectral index approaches to large values the distribution tends to the Maxwellian distri- 
' bution. In this sense, the kappa distribution is a generalization of the Maxwellian distribution. Vasyliunas 
appears to have been the first to employ the general form of the kappa distribution and to note its relation to 
the Maxwellian [1]. Kappa distribution have been used to analyze and interpret spacecraft data on the earth's 
magnetospheric plasma sheet [2] , the solar wind [3] , Jupiter [4] , and Saturn [5] . In practice it is found that many 
space plasmas can be modeled more effectively by a superposition of kappa distribution than by Maxwellians. 
In the context of both space and laboratory plasmas Hasegawa et. al. [6] showed that the equilibrium state of 
the distribution function for a plasma immersed in supcrthermal radiation resembles a Lorentzian-type distri- 
bution. The Maxwellian and kappa distributions differ substantially in the high energy tail, but the differences 

5 . become less significant as n increases. 
CO . 

Although the kappa distribution has been known for several decades, it has been mostly used for the study 

of linear waves. Under linear approximation, kappa distribution employed to study the dispersive properties 
as well as the rate of Landau damping for a number of wave modes in space plasmas such as the electron 
plasma waves, ion acoustic waves, and electromagnetic R-mode and L-mode waves and noted that the presence 
Qh of a high energy tail leads to a significant change in the damping rate of the interested wave compared to 
that of Maxwellian plasma [6,7]. However, the linear approximation can break down because of various well- 
known instability mechanisms leading to the nonlinear regime. Thus, it provides a motivation for studying 
the nonlinear effect in a plasma medium that its equilibrium state is non-Maxwellian. Note that the plasma 
equilibrium always has essential influence on the nonlinear dynamics. Nonlinear effects either directly depends 
on the velocity distribution, such as particle trapping or is affected through the average of the distribution 
function. 

In this paper, the kappa distribution function is used in the study of the ion-acoustic soliton. The ion-acoustic 
solitons in a plasma with cold ions were initially investigated by Sagdeev [8] . He assumed that the distribution of 
the electrons in the field of the wave has the equilibrium (Boltzmann) form. However, in the cases that particles 
trapping by the potential field of the wave is possible, the arrangement of particles in the phase space drastically 
changes and the closed trajectories are obtained. These particles cause a large deviation of distribution function 
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in comparison to the case when the trapped particles are not considered. Here, we follow the well-established 
kinetic model of trapping in an adiabatically varying field [9,10]. 

The present paper is devoted to a detailed calculations of the electron density from the kappa distribution 
function, including the effects of electron trapping. Since the density can not be expressed in terms of simple 
functions, therefore assuming a small potential energy (in comparison to the thermal energy), a Taylor expansion 
of the electron density is obtained. The electron density together with the fluid equations for the ions and the 
Poisson's equation form a complete set of equations that are used to study the ion-acoustic solitons. It is shown, 
as the most important result, that the deviation from Maxwellian distribution makes it possible to construct 
the localized structures which move subsonically. To our knowledge this result is quite new. Moreover, the 
soliton velocity and width decrease according to the measure of deviation. 

The lay out of the paper is as follows. In Sec. II we give the basic mathematical calculation of the electron 
density. In Sec. Ill the complete formalism and the results are presented. Section IV is devoted to the summery 
and conclusions. 

II. THE ELECTRON DENSITY 

Let us begin with the one dimensional kappa distribution function, f K (Fig. 1), for a free system as follows: 

n i r(« + i) / v 2 y K 

f K (v) = -=- — r l + -To , (1) 



v^6U 3 / 2 r(K- 1/2) V k9 2 
where k is the spectral index; the thermal speed 9 is related to the particle temperature T by 

fl=[(2«-3)//s] 1 / 2 (T/m) 1 / 2 (2) 
when k > 3/2; T is the gamma function; and f K has been normalized so that / f K dv = n - 




FIG. 1. The Maxwellian and non-Maxwellian distributions. The small rectangular shows the high energy tail associated 
with the non-Maxwellian distribution. 



We note that as k — > oo, f K — ► /m, where /m is the Maxwellian distribution function given by 

,,2 



/m = 



n 



v 



6Xp( -2# 



(3) 



where v\ = T/m. 
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We need to define an upper limit for k above which f K and /m, up to the desirable accuracy, are almost the 
same. The width of the distribution functions can be considered as a convenient choice for this purpose. We 
keep the error of the calculation of the order of one percent, and look for a solution of the following equation in 
terms of n: 



x 100=1, (4) 



where 

A fM = 2V2v T , (5) 

and 

Af K = 20 V « [exp(l//c) - 1], (6) 

are the width of Maxwellian and non-Maxwellian distribution functions, respectively. The solution of the above 
equation results in K max = 50 for the maximum value of n. 

The above equilibriums are now extended to the case when the electrons confront a well-shape potential 
energy which modifies the distribution function in the following manner: 

n i r(« + i) / v 2 -vi y K 

with the following definitions: 

n = 2 f K (v, 0) dv, (8) 
Jo 

v{ = 2U/m, (9) 

where U is just the absolute value of the potential energy. 

As it was mentioned earlier, in a collisionlcss plasma the distribution of the trapped electrons may differ 
significantly from the equilibrium distribution, which influences the properties of the wave. The particles 
density, including the trapped particles, in a Maxwellian plasma in which the fields evolve adiabatically is 
[9,10]: 

f°° 2 v 

n(v 1 ) = 2 f M (v^^dv + no—== — . (10) 
J Vl V2tt v t 

Accordingly, in the case of kappa distribution the modified number density becomes: 

"<-> - 2 [ ! " {v - n,dv + 2 ^l Jrl + -l/2r <"» 

Since the above integral can not in general be evaluated in terms of simple functions (it can be generally 
expressed in term of hypergeometric functions), we expand f K (v, vi) for small Vi(= y/2U/m) in terms of Taylor 
series and obtain 

n e . . (U\ „ fU\ 3/2 fU^ 2 

where 



nQ =l + aA- -PA- +7k I-] , (12) 



U 3^F( K -3/2) 5 / 2 r( K -3/2)' [ ' 

_l ( K +l/2)( K -l/2) 
7K_ 2 ( K -3/2) 2 ' (5) 



3 



We note that we have retained, in the expansion, the term that is proportional to U 2 . It is now possible to 
calculate the density of Eq. (11) numerically and compare it with Eq. (12)(Fig. 1). 




FIG. 2. The normalized electron density versus the normalized potential energy. 



The deviation of the expansions from the exact density is clear. The expansion is down up to U 2 for the 
following reason. Quantitatively, we define the error of the expanded density in the following form, 



error = 



exact ^-e expansion 



x 100. 



exact 

Figure 3 depicts the error for different ks versus the normalized potential energy. 
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FIG. 3. The error versus the normalized potential energy. 
It is clear that by increasing the k the error decreases. Therefore, in order to keep the error of the order of 



one percent, the U 2 term has to be kept in the density and the upper limit of the normalized potential energy 
should be restricted to 0.1. 



III. ION-ACOUSTIC SOLITON 

The ion motions are governed by the following standard equations for cold ions: 

af + T X M = °' (17) 

dvi d d , ,„ . 

-dt +Vi di Vi = -^ (l8) 

d 2 4> , . 

■^=n e -m, (19) 

where n e and rij are the electron and ion densities, vi is the ion velocity, and (f> is the ambipolar potential due 
to charge separation that have been normalized with the quantities as follows: 

Vi 

> Vi, 

6 4^<l>. (20) 
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Here e is the magnitude of the electron charge, no is the equilibrium density, c s (= y^T/rrii) is the ion-sound 
velocity, Xd[= (47re 2 no)] is the Debye length, o;^ = 47rn e 2 /m,, where m, and T are the ion mass and 
electron temperature, respectively. 

Now, assume the stationary case when all the quantities depend on £ = x — Mt where M is the velocity of 
the stationary profiles. Localized solutions are considered, with the following conditions that as £ — > oo: 

n e ,i — > 1, </> and d4>/d£ — > 0, — > 0. (21) 

Then Eq. (17) and (18) give 

20 N - 1/2 



^=^-^J • (22) 

It is worthwhile mentioning here that in order to avoid the hydrodynamic shock, the condition (f> m < M 2 /2 
has to be fulfilled. For the 4> m < 0.1, the results (Fig. 4) show that the above condition is satisfied over whole 
velocity range. Then, the ion density can be expanded for the shallow well {<f> <C 1) as follows, 

6 3 6 2 

^ = 1 + W + 2W- (23) 

Substituting the electron density from Eq. (12), ion density from Eq. (23) into Eq. (19) and integrating once 
will result to the following equation: 



1 

2 W£ 



V[<t>) - 0, (24) 



where the conditions (21) are used to determine the constant of integration. Equation (25) is the kappa 
counterpart of the Gurevich potential [9] for <f) <C 1. In fact, it is reduced to the Gurevich potential if one 
substitute the asymptotic values of a K = 1, /3 K = 4/(3v / ti : ), and j K = 1/2 in Eqs. (13), (14), and (15) when 

K — > OO. 
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FIG. 4. The ion-acoustic soliton velocity versus its maximum amplitude. 



For a soliton solution the effective potential of Eq. (25) should fulfill certain conditions. As it is seen, V((f>) 



and its first derivative vanish at (f> — and further, V(<f) has to be well-shaped between < 
is the nonzero root of V(4>) = 0. That means, a K > 1/M 2 and 



< 



where ■ 



M = 




(26) 



Figures 4 and 5 show the behavior of M 2 



versus k £ 



Mid 4> m . As it is seen from Fig. 4, solitons with larger 



amplitude move faster. Therefore, since the largest amplitude is 0.1, the upper limit of the soliton's velocity 
can be defined as follows: 



< M < 



r, - . -ir, ■ l()n. ;,. • „ 



(27) 
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FIG. 5. The ion-acoustic soliton velocity versus k. 

It is very important to note that, the deviation form Maxwcllian distribution (finite n) makes it possible 
to construct ion-acoustic solitons in the subsonic regime (Figs. 4, 5, and 6). As it was mentioned earlier, 
the velocity of ion-acoustic solitons is a monotonically increasing function of its maximum amplitude. In the 
amplitude interval over which the above theory is based (0 < 4> m < 0.1), as it is seen from Figs. 4 and 5, 
that solitons, resulting from small n (k < 6) are subsonic. Solitons associated with the larger k, according to 
their maximum amplitude move sub or supersonically. Figure 6 shows the boundary of subsonic and supersonic 
regimes (M = 1). 
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FIG. 6. The sub and supersonic regimes. 

For larger k the amplitude interval over which the solitons are subsonic gets narrower. However, the existence 
of subsonic soliton is particular result of non- Maxwcllian plasmas. Solitons move supersonically all over the 
amplitude interval just when n approaches to infinity, or in the other words, when the plasma electrons become 
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Maxwellian. 

We now discuss the width of the soliton. The width of the soliton is defined as 

ra max /e' da 

A = / ; (28) 

where e' in the upper limit is the Neperian number and the maximum amplitude has been assumed to be at 

e = o. 
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FIG. 7. The sub and supersonic regimes. 

Figures 7 and 8 exhibit the dependence of the width A on the maximum amplitude (f) m and k. Apparently, for 
any <j> m the width of the solitons resulting from non-Maxwellian deviation is narrower according to the measure 
of its deviation from Maxwellian equilibrium. 
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FIG. 8. The sub and supersonic regimes. 
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IV. SUMMARY AND CONCLUSIONS 



In this paper, the nonlinear theory of plasma waves was extended to the plasmas that their equilibrium state 
are specified by the non-Maxwellian (here kappa) distribution. We believe that the extension is very important 
since most of the space and some of laboratory plasmas are not in the Maxwellian equilibrium. Although the 
linear theory of this issue has been known for the decades but, to our knowledge, this is the first attempt in 
opening the gate to the nonlinear world of plasma waves in the non-Maxwellian background. As an example the 
ion-acoustic solitons are studied in this framework and it was shown that as the most important result regarding 
to the non-Maxwillian equilibrium is the possibility for solitons to move subsonically. The solitons velocity and 
width arc monotonically increasing function of spectral index (k) and approach to the their Maxwellian values 
as k — > oo. 
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